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Abstract 



The formation and late time decays of Q-balls are generic consequences of the 
Affleck-Dine (AD) baryogenesis. A substantial amount of the lightest super- 
symmetry (SUSY) particles (LSPs) are produced non-thermally as the decay 
products of these Q-balls. This requires a significantly large annihilation cross 
section of the LSP so as not to overdose the universe, which predicts a higgsino- 
or wino-like LSP instead of the standard bino LSP. We have reexamined the 
AD baryogenesis with special attention to the late-time decays of the Q-balls, 
and then specified the parameter regions where the LSPs produced by the Q- 
ball decays result in a cosmologically interesting mass density of dark matter 
by adopting several SUSY breaking models. This reveals new cosmologically 
interesting parameter regions, which have not attracted much attention so far. 
We have also investigated the prospects of direct and indirect detection of 
these dark matter candidates, and found that there is an intriguing possibility 
to detect them in various next generation dark matter searches. 
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I. INTRODUCTION 



Both the origin and nature of dark matter and the production mechanism of the baryon 
(matter-antimatter) asymmetry in the present universe are among the most fundamental 
puzzles in cosmology and in particle physics. Of the many candidates for dark matter, 
the lightest supersymmetric particle (LSP) has been attracting wide interest as an ideal 
candidate, which is inherent in the supersymmetric (SUSY) standard model, and it is exactly 
stable under the i?-parity conservation. 

On the other hand, once we introduce SUSY in the standard model, there appear a 
number of flat directions in the scalar potential carrying baryon (B) and/or lepton (L) 
number. Affleck and Dine (AD) proposed an interesting mechanism to generate effectively 
a baryon asymmetry by using these flat directions. Because of the flatness of the potential, 
it is quite plausible that a linear combination of squark and/or slepton fields has a large 
expectation value along a flat direction during the inflation in the early universe. After the 
end of inflation, this flat direction field <fi starts a coherent oscillation. At this stage, due to 
baryon (and/or lepton) number- violating operators in the scalar potential, the (f) field gets a 
non-zero motion along the phase direction as well, which means nothing but baryon (and/or 
lepton) number generation since the (f) carries a baryon (and/or lepton) number. 

Recently, we have pointed out that the higgsino- or wino-like neutralino naturally becomes 
the dominant component of the dark matter, if either of them is the LSP and if the AD 
mechanism is responsible for the generation of the observed baryon asymmetry in the present 
universe The LSP is produced non-thermally by the late-time decay of the Q-balls, which 
are generally formed in the AD mechanism 0-0. Then, the large pair annihilation cross 
section of the higgsino- or wino-like LSP leads to the desired mass density of dark matter, 
instead of overdosing the universe, as in the case of the standard bino LSP. It should be 
noticed that the typical decay temperature of the Q-ball is between 1 MeV and (a few) GeV, 
i.e. below the freeze-out temperature of the LSP and before the big-bang nucleosynthesis, 
which is just the desired range for the present scenario. In this scenario, it is very interesting 
that both the dark matter and the baryon asymmetry are produced by a single source, the 
Q-ball decay, which originates from the AD condensation. Furthermore, the relatively strong 
interactions of these LSPs significantly enlarge the possibility of detecting these dark matter 



1 We do not consider the gauge-mediated SUSY breaking models ||, where the Q-ball is generally 
stable |||. 
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candidates in both direct and indirect searches. 

In this paper, we give more detailed analyses of this scenario, "higgsino or wino dark mat- 
ter from Q-ball decay" , adopting several SUSY breaking models, which include the minimal 
supergravity (mSUGRA) scenario, the anomaly-mediated SUSY breaking (AMSB) model |§ 
with additional universal scalar mass, and the no-scale model with non-universal gaugino 
masses 0. In each model, we will show parameter regions in which the LSP from the Q-ball 
decay can be the dominant component of the dark matter. This predicts new cosmologi- 
cally interesting parameter regions, which have not attracted much attention; the wino dark 
matter is realized in wide parameter regions in the AMSB models and the no-scale models, 
where the non-thermal wino production via decays of Q-balls naturally explains the mass 
density of the dark matter. Even in the mSUGRA scenario, higgsino dark matter is realized 
in a relatively wide region by virtue of the "focus point" behaviour of m 2 H |K| , in which the 
mass density of dark matter is explained by the non-thermally produced higgsino. 

We also investigate the direct and indirect detection of the neutralino dark matter in these 
regions. We calculate the proton-neutralino scalar cross section for the direct detection. As 
for the indirect detection, we adopt the XX ~* 77 annihilation channel, which produces 
monoenergetic 7-ray lines. We estimate the neutralino annihilation rate into the 27 final 
state using the full one-loop calculations presented in Ref. Hll| , and discuss the possibility 
of detection in the next generation of air Cherenkov telescopes for several models of density 
profile in the halo. Actually, the large higgsino or wino content of the LSP strongly enhances 
the direct and indirect detection rates, via the Higgs exchange diagrams and the chargino-W 
boson loop diagrams, respectively. As we will see, if the higgsino or wino is really a significant 
component of dark matter, there exists an intriguing possibility to detect these non-thermal 
dark matter candidates in various next-generation detectors. 

Note that, if the higgsino- or wino-like dark matter is indeed detected at future exper- 
iments, it suggests the existence of non-thermal sources for these LSPs, which predicts a 
highly non-trivial history of the universe. This might discriminate the origin of the baryon 
asymmetry between various scenarios. 

The rest of this paper is organized as follows: first, we reexamine the AD baryogenesis 
scenario, giving particular attention to the late-time decay of Q-balls in Section |TJ. This 
includes an interesting new twist in the AD baryogenesis using non-renormalizable operators 
in the Kahler potential to rotate the AD field. This scenario is now able to explain the 
observed baryon asymmetry without any assumptions of additional large entropy production 
to dilute the resultant baryon asymmetry, such as decays of heavy moduli fields and/or a 



3 



thermal inflation. In Section |HI| , we review the properties of the Q-ball in the minimal SUSY 
standard model (MSSM). We discuss the non-thermal production of the LSP dark matter 
from the Q-ball decay in Section IV. Our main results are given in Section We show 



in several SUSY breaking models the parameter spaces where the higgsino- and wino-like 
LSPs from the Q-ball decay can be the dominant component of the dark matter, and we 
investigate their direct and indirect detection in detail. We also comment on the possibility 
that the non-thermally produced neutralino via the late-time decays of Q-balls form the 
warm dark matter, which erases the cuspy profile of the matter density of the halo that 
might be inconsistent with the observations. We present concluding remarks and discussions 



in Section VI. 



II. AFFLECK DINE MECHANISM 

In this section we briefly review the Affleck-Dine mechanism HIJ12fl in the MSSM. Here- 
after, we assume i?-parity conservation and the following superpotential at renormalizable 
level as usual: 

W = yu QUH u + y D QDH d + y E LEH d + fiH u H d . (1) 

Here, Q, U, D, L, E, H u and H d denote superfields of left-handed quark doublets, right- 
handed up-type and down-type quarks, left-handed lepton doublets, right-handed charged 
leptons and Higgs doublets, respectively; we omit the family indices for simplicity. In the 
scalar potential, there are a number of flat directions along which the F-term potential coming 
from the superpotential in Eq. (|T]) as well as the D-term potential vanish : 

LH U , H u H d , 

UDD , LLE , QDL , 

QQQL , UUDE , QUQD , QULE , • • • , (2) 

where the ellipsis denotes the flat directions consisting of linear combinations with more 
scalar fields. In the following, we will parametrize a flat direction by a complex scalar field 
0. 

A. Affleck Dine mechanism with non-renormalizable operators in the superpotential 

In the supersymmetric limit, the flat directions are lifted only by non-renormalizable 
operators in the scalar potential. Let us first assume that all the non-renormalizable operators 
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consistent with the MSSM gauge symmetry and the i?-parity exist in the superpotential. 
Most of the flat directions are then lifted by the scalar potential coming from the following 
dimension-4 superpotentials: 

W = ^LH U LH U , ^QQQL , ^UUDE , (3) 

where M denotes effective scales at which these non-renormalizable operators are induced. 
Here and hereafter, we discard superpotentials that conserve the baryon and lepton numbers 
(e.g. W oc QUQD, QULE), since they cannot generate baryon asymmetry. We have included 
superpotentials conserving B — L, since the Q-ball can decay after the electroweak phase 
transition and the baryon asymmetry is not washed out by the sphaleron effect |14] . 



Some of remaining flat directions are lifted by the following dimension-6 superpotentials: 

W = -^-UDDUDD ,-^-LLELLE . (4) 

M 3 M 3 w 

All the other flat directions are lifted either by B- and L- conserving superpotentials or by 
superpotentials of the form W = (X/M n ~ 3 )if)(j) n ~ 1 [I3"|. In the latter case, one can show that 



t/> becomes zero when the <fi field develops a large vacuum expectation value. Therefore, this 
type of superpotential cannot provide a non-zero A-term that is indispensable to create a 
net baryon/lepton asymmetry as will be discussed soon. Thus, in the following, we write the 
superpotential that lifts a flat direction <fi as 

with n = 4 or 6. Note that we have normalized the effective scale M including the coupling 
A, and hence M can be larger than the cut-off scale (e.g. Planck scale). 

In the case of the leptonic flat direction (e.g. LH U , LLE), lepton asymmetry is first 
generated by the AD mechanism. Then a part of the produced lepton asymmetry is con- 
verted [fn| into baryon asymmetry through the sphaleron effect. In particular, leptogenesis 
via the LH U flat direction |]I6| has attracted much attention [fT7-22||, since in this case the 
baryon asymmetry in the present universe is directly related to the neutrino mass. This al- 
lows us to obtain definite predictions on the rate of the neutrinoless double beta decay [|I~9| , [2i| . 
In this paper, however, we do not consider the LH U flat direction since the resultant Q-ball 
along this direction is very small and evaporates well above the weak scale P,|25[|. As for 
the other leptonic flat directions, large Q-balls can be formed. However, the Q-balls must 
decay or evaporate well before the sphaleron effect terminates (T ~ 100 GeV), which requires 
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some mechanisms to make Q-balls small enough [ 24], 25j] . In this case, the LSPs produced by 
the decays of Q-balls are thermalized and the bino-like LSP is the almost unique candidate 
for neutralino dark matter. We do not treat such cases in the following discussions, and 
concentrate on the flat directions carrying non-zero baryon number. 

Let us now discuss the AD mechanism with non-renormalizable superpotential. The 
relevant scalar potential for <j) is given by 

v(0) = K- c ^ 2 )i0i 2 

+^(^ + ^-) + ^l0| 2 "- 2 - (6) 

Here, the potential terms, which are proportional to the soft mass squared and the grav- 
itino mass 7713/2, come from the SUSY breaking at the true vacuum. We will concentrate on 
gravity-mediated and gaugino- mediated SUSY breaking models, and take ~ rn 3 / 2 \a m \ ~ 
1 TeV, hereafter. (We will discuss later the case of anomaly-mediated SUSY breaking mod- 
els, where m 3 / 2 ^> m^.) The terms depending on the Hubble parameter H denote the effect 
of SUSY breaking caused by the finite energy density of the inflaton [12] : ch is a real con- 
stant of order unity, which depends on the couplings between the inflaton and the <fi field 
in the Kahler potential. Hereafter, we take cjj — 1 (> 0), which is crucial to let </> have a 
large expectation value during inflation. There might also exist a Hubble-induced A-term 
potential, which has the same form as the second term in Eq. (|6|), with H instead of m3/2.Q 
Even in the presence of this term, the following discussions are not altered^] so we assume 
its absence in the present work, for simplicity. The last term is the F-term potential coming 
from the superpotential in Eq. (H). Here, we assume the absence of thermal effects, which 
will be justified later. 

The estimation of the baryon asymmetry is rather straightforward. The baryon number 
density is related to the AD field as 



2 This requires the existence of a three-point coupling of the AD field to the inflaton in the Kahler 
potential, 5K D I<jy(j)/Mpi + h.c, where / denotes the inflaton superfield. This coupling leads to 
a relatively high reheating temperature of the inflation, Tr ~ m/(m//M p j) 1 / 2 , where mi is the 
mass of the inflaton. Therefore, as we will see, the absence of this three-point coupling is also 
desirable for obtaining the low reheating temperature that explains the right amount of baryon 
asymmetry. Notice that ch — 1 does not need this three-point coupling, which is consistent with 
the low reheating temperatures. 

3 See the related discussions in Refs. [12,17,19|. 
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n B = f3i{<p*(f) - 4>*4>) , 



(7) 



where j3 is the corresponding baryon charge of the AD field, which is at most 1/3. The 
equation of motion for the AD field in the expanding universe is given by 

^ + 3^0+^ = 0. (8) 

Then, together with Eq. (0), the equation of motion for the baryon number density is written 
as follows: 

fdV \ 

h B + 3Hn B = 2(3 Imi—cj)] 

= 2(3 ^ Im (a m r) , (9) 

where we take 777.3/2 as real by adjusting the phase of a m . One can see that the production 
rate of the baryon number is proportional to the A-term. By integrating this equation, we 
obtain the baryon number at the cosmic time t as 

R 3 n B ] (t) = 2(3 f R 3 \<p\ n sine dt , (10) 

where R is the scale factor of the expanding universe, and 9 = arg(a m ) + 77 arg(0). 

During inflation, the negative Hubble mass term — c#if 2 |0| 2 causes an instability of the 
flat direction field <fi around the origin, and the AD field acquires a large expectation value: 

|0| ~ (ff;M"- 3 ) 1/( " 2) , (11) 

where Hi denotes the Hubble parameter during inflation. This is the balance point between 
the negative Hubble mass term -c H H 2 \(f)\ 2 and the F-term potential \^\ 2n ~ 2 j M 2 ^ . Note 
that the curvature along the phase direction m 2 h is much smaller than iff at this point, since 

m lh * ^Sl^r 2 ^ m ^Hi « Hj . (12) 

Therefore, the initial phase of the AD field fixed during inflation does not generally coincide 
with the minimum of the A-term potential in Eq. (|6|), and hence we naturally expect that 
sin^ = 0(1). 

After the end of inflation, the AD field slowly rolls down toward the origin following the 
gradual decrease of the Hubble parameter H as \<f>(t)\ — (H(t)M n ~ 3 ) 1 ^ n ~ 2 . At this slow rolling 
regime, the right-hand side of Eq. ( |10"D increases as oc t 2 ~ 2 ^ n ~ 2 \ Here, we have assumed the 
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matter-dominated universe, which is true as long asT fl <2x 10 10 GeV(m /lO 3 GeV) 1 . 
When H(t osc ) ~ m^, the soft mass term of the AD field eventually dominates the negative 
Hubble mass term, and causes the coherent oscillation of the AD field around the origin. 
After this time, the amplitude of the AD field rapidly decreases as \<f>\ oc t~ x , and then the 
production of the baryon number terminates at the time H osc ~ m^. 

Using the above arguments and Eq. ([II]), we obtain the baryon number density at the 
time t = t osc : 

2(n — 2) / „\2/(n-2) 

^B(tosc) = 3 n _ 3 } , /35 eS \a m \m 3/2 (ff osc M"- 3 ) , (13) 

where 6 e g = sin#(= 0(1)). Then, after completion of the reheating process of the inflation, 
this leads to the following baryon asymmetry: 

n B 1 T R 



s 4^ 



n(t osc ) (14) 



U ~ 2 at ,„ | m 3/ 2 T fl / 3 N2/(n-2) 



where Tr is the reheating temperature of the inflation, s is the entropy density of the universe 
and M p i = 2.4 x 10 18 GeV is the reduced Planck scale. By using H osc ~ m^, it is given by 

(16) 

for n = 4, and 

=2 - 5 x 10-" x /^M (Tm) fi^XV' 2 (il) 3 ' 2 , (17) 

for n = 6. These quantities remain constant unless there are additional entropy productions. 

These predictions on the baryon asymmetry are not altered by the presence of the Q-ball 
formation, as long as the produced Q-balls are unstable and decay before the big-bang nu- 



cleosynthesis (BBN). The details related to the Q-balls will be discussed in Section |IJ. The 
only point remaining to be checked is the absence of the thermal effects. In the remainder of 
this section, we discuss these in the cases of the n = 4 and n = 6 superpotentials. 

1) n = 4 

In the case of the n = 4 superpotential, it was pointed out that the thermal effects 
cause the early oscillations of the AD field H osc ^> m^, which significantly suppress the 
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resultant baryon asymmetry p6 . Analytic and more systematic investigations were done 
for the case of the leptogenesis via the LH U flat direction, which revealed an interesting 



property, the "reheating temperature independence" of the baryon asymmetry ||19|| . For the 
other flat directions, however, there are some subtleties. The dimension-5 operators coming 
from the superpotential oc QQQL, UUDE are responsible for the proton decay [27j. Null 



observations of the proton decay then lead to strong constraints on the effective cut-off 
scale, 10 25 GeV, at least for the most relevant operators |2^j29|] . Following the analyses 
in Ref. fT9j, one can show that there is no thermal effect as long as Tr^S 10 6 GeV when 
M ^ 10 23 GeV, and hence there is no need to worry about the thermal effects in these 
scenarios. 

As we will see, such a large scale of M leads to the formation of very large Q-balls with 
the decay temperature T d ^ 1 GeV. Since the decays of the Q-balls occur well after the 
sphaleron effect terminates, the baryon asymmetry is not washed out, although the relevant 
flat directions conserve the B — L symmetry. 



2) n = 6 

In the case of the n = 6 superpotential, it is much easier to avoid the thermal effects 
because the potential is much flatter than in the n = 4 case. For the appearance of the 
thermal mass terms, the fields coupled with the <fi field must reside in the thermal bath: 
/| 0| < T, where / is the Yukawa or gauge coupling constant of the AD field. Then, the 
sufficient condition to avoid the thermal mass terms is written as 



T R < 10 8 GeV 



/ 



M 



3/2 



v io- 5 ; v m pi/ 

Another thermal effect we have to check comes from the thermal logarithmic potential |30 



18) 



5V D aT 4 log 



J>2 



(19) 



where a is a constant given by the fourth power of gauge and/or Yukawa coupling constants. 
To avoid the early oscillations due to this potential, we need the following condition: 



4 Another possibility is the absence of these operators. In such a case, most of the flat directions 
are expected to be lifted by the n = 6 superpotential, and two other operators are operative for AD 
baryogenesis besides (UDD) 2 , such as UDDQDL, LLEUDD. 
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Tr < JL M (_£U 1/4 = 7 . x 10 r GeV tlgS 1/2 < MS m ( -t_f , (20) 

which leads to the comparable constraint on the reheating temperature as Eq. (|T8|). As 
a result, the early oscillations due to the thermal effects can be easily avoided as long as 
Tr J$ 10 8 GeV in AD baryo/leptogenesis with n = 6 superpotential. 



B. Affleck Dine mechanism without superpotential 



The other interesting possibility is that non-renormalizable operators in the superpotential 
are forbidden by some chiral symmetries, such as i?-symmetry. In this case, the initial 
amplitude of the AD field is determined by the negative Hubble mass term — cnH 2 \(f)\ 2 and 
the non-renormalizable operators in the Kahler potential and is expected to be |0|o ~ M p \. 

In this section, as an example, we consider the following terms JE], which are consistent 
with i?-symmetry: 

r (7^7 - 7^ 7 \ 

5C = j d 4 9 ^—QQWE^ + \ 2 —QUW^L + h.c.J , (21) 

where Aj's are coupling constants and Z is the superfield with the non-vanishing F-term 
responsible for the SUSY breaking in the true vacuum. The effective potential of the AD 
field induced by these terms and SUSY breaking effects is then given by 

2 

^ = K-c^ 2 )|0r + ^(A0 4 + h.c.) + ... , (22) 

pi 

where A is a coupling constant and the ellipsis denotes the higher order terms coming from the 
Kahler potential. (Above the Planck scale the scalar potential is expected to be exponentially 
lifted by the supergravity effect.) 

The evolution of the AD field is much simpler than the case with non-renormalizable 
superpotential. During inflation, the AD field is fixed at a very large scale |0|o ~ M v \ because 
of the negative Hubble mass term —chH 2 |0| 2 . The AD field just stays there until H = H osc ~ 
77i0 because of a large damping effect, which appears in the second term in Eq. (^. The baryon 
number production completes as soon as the AD field starts the coherent oscillation around 
the origin at H = H osc . In this scenario, the AD field is completely decoupled from thermal 
backgrounds, and the early oscillation of the AD field (H osc ^> m^) thus does not occur. 

A crucial difference appearing in this scenario is the existence of a large entropy production 
due to decays of the AD field. At the reheating process of the inflation, the energy density 
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of the inflaton is converted into the radiation. After that, the energy density of the radiation 
reduces as p ra d oc R~ A . On the other hand, the AD field starts the coherent oscillation 
at H = H osc with the initial amplitude \<p\ = \<p\ , and its energy density only reduces as 
P0 oc R~ 3 . Then the AD field begins to dominate the energy density of the universe at some 
time, and substantial entropy is produced through decays of the AD field. 

The AD field dominates the energy density of the universe before its decay if the following 
condition is satisfied: 

T R > 3 T d , (23) 

where Td is the decay temperature of the (f> field, which corresponds to the decay temperature 
of the Q-ball and much lower than the weak scale, Td <C m w . In Eq. (|23]), we have assumed 
that the coherent oscillation of the AD field starts before the completion of the reheating 
process of the inflation. If this is not the case, i.e. if T R ^ 2 x 10 10 GeV (m^/l TeV) 1 / 2 , the 
condition on Td becomes weaker than this case. 

After the decays of the AD field, the resultant baryon asymmetry is given by the following 
simple form: 

nB_ _ Pi f nA ( niA _ 3 / n B \ 
s s \ / \ ) 4 \ ) 

where denotes the total number density of <fr and anti-0 field. The ratio of baryon to <fi 
number is easily estimated by an argument similar to the one given in the previous section 

as 

2 

I S eS , (25) 




J \Mp\ / 

where 5 e s = sin(arg(A) + 4arg(0 o ))- Then a reasonable set of parameters 

can naturally explain the observed baryon asymmetry. 

Before the discovery of the formation of large Q-balls, decays of the AD field have been 
considered to occur at Td ~ m^. If this is the case, the resultant baryon asymmetry is washed 
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out by the sphaleron effect, since the Kahler potential adopted in the present work conserves 
the B — L symmetry. If we adopt other interactions that violate this symmetry,^ the resultant 
ratio of 0-number to entropy is given by n^/s ~ 0(1), which requires a very small ratio of 
baryon to 0-number Ub/u^ ~ 10~ 10 . Such a small ratio results in a very unnatural fine tuning 
on A and/or 5 e g, and hence additional large entropy productions have been considered to be 
necessary after the decays of the AD field, such as a thermal inflation or decays of heavy 
moduli fields. 

Interestingly enough, the large initial amplitude of the AD field results in the formation of 
very large Q-balls with low decay temperature Td = O(10 MeV).f] This significantly enhances 
the entropy production via decays of the AD fields and naturally explains the observed baryon 
asymmetry. Higgsino or wino LSP is the necessary condition for this scenario to work, so as 
not to overdose the universe and explain the right amount of dark matter. 

A quite beautiful point in this scenario is that the resultant baryon asymmetry and dark 
matter density are solely determined by the sector related to the AD field and annihilation 
cross section of the LSP. Any other process cannot affect the final result as long as the 
AD fields (stored in the Q-balls) dominate the energy density of the universe at their decay 
time. Especially, the final baryon asymmetry is completely independent of the reheating 
temperature of the inflation. 



In addition, the so-called "cosmological gravitino problem" [|3T] is also solved. The dilution 
factor, which is the ratio of the entropy density before and after the decays of the AD field, 
is given by 

i 5*. = 3^L^y (27) 

A s after T ini \J0|o/ 

where 

/ ( 90 \ 1/4 " 

T ini = Min T R , y/m^M^ — , (28) 

with being the effective degrees of freedom at temperature T. It can be seen from these 
relations that even if the reheating temperature of the inflation is much higher than ~ 
10 12 GeV, the produced gravitinos are sufficiently diluted for there not to be any gravitino 
problem. 



5 Such interactions consistent with i?-symmetry only appear at much higher orders. 
6 See Sec. EH. 
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C. Affleck Dine mechanism in anomaly- mediated SUSY breaking models 



In anomaly-mediated SUSY breaking (AMSB) models, SUSY is broken in a hidden sector, 
and it is transmitted to the MSSM sector dominantly via the super- Weyl anomaly. The 
MSSM gauginos and matter fields obtain soft masses, which are one-loop-suppressed relative 
to the gravitino mass. Therefore, the order of the gravitino mass is estimated as m 3 / 2 ~ 
m soit /a > TOgoft. 

A crucial point is that the AD mechanism always uses the non-renormalizable operators 
to violate baryon and/or lepton number, which also violate the super- Weyl symmetry at tree 
level. Hence, the SUSY breaking effects induced in these operators are of the order of the 
gravitino mass and not loop-suppressed, which generates a global minimum for the AD field 
displaced from the origin 



First consider the AD mechanism using non-renormalizable operators in the superpoten- 
tial. In this case, the scalar potential of the AD field is given by Eq. (Eft, but with much 
larger gravitino mass. During the inflation, the negative Hubble mass term drives the AD 
field far from the origin and fixes it at ~ (f/j-M™ -3 ) 1 /™ -2 . After the inflation ends, the 
amplitude of the AD field gradually decreases as |0| ~ (ifM™ -3 ) 1 /™ -2 . At this stage, the 
curvature along the phase direction m^ h is given by m^ h ~ m 3 / 2 iJ. Thus, after H^m 3 / 2 , 
this large curvature forces the AD field to settle down on the bottom of one of the valleys of 
the A-term potential, which are located at 

arg(a m ) + n arg(0) = ir, mod n . (29) 

Consequently, the AD field is inevitably trapped at the global minimum on the way toward 
the origin along the bottom of the valley: 

/ Q \l/(n-2) 

|0| min ~ (m 3/2 M"- 3 ) . (30) 

This results in an SU(3)c and/or U(1)em breaking universe. To avoid this disaster, we need 
some mechanisms to stop the AD field during the inflation below "the top of the hill" in the 
scalar potential, which is located at 



lo 



< 



m 



2 \ V(n-2) 

rn— 3 



M n ~ A . (31) 



The same difficulty appears in the AD mechanism making use of the non-renormalizable 
operators in the Kahler potential. In this case, the potential is written as in Eq. (0), with 
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larger gravitino mass m 3 / 2 3> m^, and the global minimum is located near the Planck scale. 
In order to avoid a colour and/or U(1)em breaking universe, we need to impose the following 
condition: 

l4><Mun = Mpr^. (32) 
m 3 / 2 



One natural solution is to gauge the U(1)b-l symmetry |24]].fl If the B — L symme- 
try is gauged, flat directions along which the AD fields cannot cancel the U(1)b-l D-term 
contributions within themselves can be lifted at the B — L breaking scale v by the D-term 
potential. Then, the amplitude of the AD field can be easily suppressed by adjusting the 
B — L breaking scale as v ^ \4>\bm- Even if we suppress the initial amplitude of the AD field, 
there remains plenty of room to explain the observed baryon asymmetry and to form large 
Q-balls, which allows us to obtain the right amount of wino dark matter in AMSB models. 

III. Q-BALL 

After the flat direction field <fi starts its coherent oscillation, a kind of non-topological 
soliton "Q-ball" |33| is formed because of spatial instabilities of the field [0-0]. In this 



section we briefly review the formation and the decay of the Q-ball. 

A. Size of the Q-ball 

The relevant scalar potential of at the time of Q-ball formation is given by 

V{<P)=ml(l + K\og(^])\<P\\ (33) 



where Mq is the renormalization scale at which the soft mass is defined, and the 
i^logd^] 2 ) term represents the one-loop correction. This mainly comes from the gaugino 
loops and K is estimated in the range from —0.01 to —0.1 P,[7|,p3| . Because of the potential 



7 In Ref . , we make use of the gauged B — L symmetry to make the Q-ball small enough for the 
produced bino LSPs not to overdose the universe. In the present work, we concentrate on the case 
of the higgsino or wino LSP, so that there is no need for such an adjustment on the B — L breaking 
scale. 
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flatter than 2 , spatial instabilities of the homogeneous <fi field are caused after it starts the 
coherent oscillation ||. In the momentum space, the instability band is given by |34 



< |r < 3mJ 



K\ 



(34) 



where k is the comoving momentum of the fluctuations of the (f) field. The best amplified 
mode is given by the centre of the band: {k 2 / R 2 ) max — (3/2)mi\K\. This corresponds to 
the radius of the Q-ball, which is estimated analytically using the Gaussian profile of the 
Q-ball 0]: 



R 



Q 




For the best amplified mode, the perturbations 
following equation: 



(35) 

of the <j) field grow according to the 





Scj) 




exp 






osc 




5<j) 


exp 







osc 



dt 



3Tn<f,\K\ 



2H 



(36) 



where we have used H~ l = (3/2)t, assuming that the Q-ball formation occurs before the 
reheating process of the inflation is completed.^ Therefore, the fluctuation becomes non- 
linear when the Hubble parameter becomes 

_ m<t,\K\ 



H = H n 



2a 



(37) 



where a = ln(|0| osc /|50| osc ) and 



(2itRq) 1 for the best amplified mode For the 



parameters we adopt in the present analyses, its typical value is given by a ~ 30-40. At this 
time (H = H non ), the baryon number density of the AD condensation becomes 

2 



n B (t r 



n B (t c 



x 



(3c 



8 In the AD mechanism without non-renormalizable superpotentials, the reheating temperature 
can be as high as Tr ^2x 10 10 GeV(m^/l TeV) 1//2 and the Q-ball formation might occur after the 
reheating phase. (See Sec. [Q B| .) In this case, H^ 1 = 2t, and Eq. (38) is modified as ns(t non ) — 
"_b(^osc) x (H non / H osc ) 3 / 2 . In such a case, the charge of the Q-ball becomes larger by a factor of 
O(10). 
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Thus, from Eqs. (0), fl37|) and ( p8|) and if osc — ^</>, the typical charge of a single Q-ball is 
estimated to 

4 

<2 ^ -vt^Rq x n B {t non ) 

2yft*\K\W 1 . , 
" "I ^^J nB(W) 

^.(MV (39) 

where 

Here, we have used n^(t osc ) ~ m^l^l^. This estimation of the Q-ball charge is roughly- 
consistent with the result obtained by numerical lattice simulations |34| , |36| . 

Q^J^) 2 xf e [° r e | £c ' (41) 
\ j [ e c for e^e c , 

where 

P n = 6x 1(T 3 , e c ~ 0.01 . (42) 

Notice that for e ^5 e c the Q-ball charge becomes independent of the initial charge density 
e. This is because negative charge Q-balls are created in this region pi" |. In the following 



discussions, we shall adopt the following equation in estimating the Q-ball charge: 

( e 



Q = Qr 



-for e c < e < (3 , 

(43) 

■ - for e < e r , 



I /3 



where 

g max ~3xlO- 3 f^V/?. (44) 

Let us now estimate the charge of the Q-ball in each of the cases discussed in the previous 
sections. In the case of AD baryogenesis with non-renormalizable superpotential (Sec. [II A|) , 
the Q-ball charge is estimated by using Eqs. (|13|) and (f43|), resulting in 
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2(n - 2) 
3(n - 3) 



P5 cS \a r 



rrir, 



(45) 



Q 



3(n-3)' 
3xl0 2O x/35 eff |a m |(^ 



1 TeV 



2(n-3)/(n-2) 



M 



10 26 GeV 



3 x 10 2U x /?<5 cff | a r , 



M 



m 



1 TeV \ 3/2 
m J VMp^ 



3/2 



for n = 4 , 
for n = 6 . 



(46) 



As mentioned in Section PB| , the Q-ball charge becomes larger in the AD mechanism 
without non-renormalizable superpotential, since the initial amplitude of the 0's oscillation 
is large. Recall that a relatively smaller value of e = ns/n^ < e c is preferred in order to 
generate the correct amount of baryon asymmetry [see Eqs. fl25|) and ( p6|) 1. Thus, from 
Eq. fl43|), the Q-ball charge is estimated to be 



Q ~ 10 



2^ 



| \ 2 (I TeV\ 2 



M, 



pi, 



(47) 



which results in Q ~ 10 24 -10 26 for |0| o ~ (0.1-1) x M pl . 

Finally, in the case of the AMSB model, the initial amplitude of the AD field should 
be suppressed to avoid the unwanted trapping of the AD field at the global minimum, as 
discussed in Section [11 (J| . From Eqs. (|31~D, ( p2|) and 

Q < 10- 3 



, the Q-ball charge is estimated to 



T) — 4 1 

18 _ ^30 TeV\ ( M 



3 x 10 18 x 



2 x 10 19 x 



m 



3/2 



/ V10 26 GeV 

1/2 



30 TeVV / 1 TeV 



M 



3/2 



M 



pi, 



"T.3/2 / V m 4> 

for the case with non-renormalizable superpotentials, and 

Q Z 10- 3 



for n = 4 . 
for n = 6 . 



(4? 



M 



pi 



, m 3/2, 



6 x 10 24 



30 TeV 

, ™3/2 , 



(49) 



for the case without non-renormalizable superpotentials. 
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B. Decay of the Q-ball 



Before discussing the decay of the Q-ball, we comment on its evaporation p7| , |38| . Because 
of the large expectation value of the scalar field inside the Q-ball, most of it is decoupled from 
the surrounding thermal bath. However, a thin outer region of the Q-ball is thermalized, since 
particles in the thermal plasma can penetrate into this region, and hence a partial evaporation 



of the Q-ball charge occurs. As stressed in Ref. ||38|| , the evaporation of the charge from the 
Q-ball surface is suppressed by diffusion effects, and the evaporation occurs most effectively 
at T ~ 777,0. The total amount of the evaporated charge is estimated to 

1/2 

\K\ 




for T R ^ m^,, and 



' 0.03 \ 1/2 (I TeV\ 3 / T fi 



for Tji^m^. (Here, we have used T ~ (T^M pl if) 1//4 before the reheating process of the 
inflation is completed.) Therefore, as long as the initial charge of the Q-ball is larger than 
O(10 18 ), most part of the Q-ball charge survives the evaporation. Hereafter, we assume 
Q ^ O(10 18 ) GeV (which is naturally realized in the present scenarios, as shown in the pre- 
vious sections) and discuss the emission of the remaining charge by the decay of the Q-ball. 
The remaining charge of the Q-ball is emitted through its decay into light fermions. The 



decay rate is estimated as [39 



r --^£^, (52) 



dQ uj 3 A 
~dl ~ 192vr 2 

where A = ^ttRq is the surface area of the Q-ball and u ~ m^. This upper bound is likely to 
be saturated for 0(0) 3> m^, where 0(0) is the field value of the AD field at the centre of the 
Q-ball |^J.[] Therefore, the lifetime of a single Q-ball with an initial charge Qi is estimated 



to be 

Tj = 

r 



r d = ^>lxl0- 7 secx [g±\[LJZL\l™\ , (53) 



Q 




3 If kinematically allowed, there are also decay channels into lighter scalars through 3-point cou- 



plings. Even if this is the case, the decay rate is at most comparable to that in Eq. (52). 
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or equivalently, the decay temperature Td of the Q-ball is given by 

Therefore, the Q-ball decays at ~ 1 MeV-(a few) GeV for the Q-ball charge estimated 
in the previous sections. It is quite interesting to observe that the Q-ball decay occurs just 
after the freeze-out of the LSP (7/ ~ m x /20) and just before the beginning of the big-bang 
nucleosynthesis (BBN) (T ~ 1 MeV) . Thus, the Q-ball decay can naturally provide both the 
baryon asymmetry required from the BBN and the non-thermal source of the LSP. It should 
also be noticed that the decay occurs after the electroweak phase transition (T ~ 100 GeV). 
Thus, the produced baryon asymmetry is not washed out by the sphaleron effect, even if the 
B — L is conserved in the AD mechanism, as in the n = 4 cases and the W = cases. 



IV. NON-THERMAL DARK MATTER FROM THE Q-BALL DECAY 

We now turn to the discuss of the non-thermal production of the LSP from the Q-ball 
decay. The relation between the Q-ball number density and baryon number density is given 
by 



n B = Qi(n^-n Q ) 



1 for e c < e < (3 , 



I _£ for e < 6c) 
e 

where Qi is the absolute value of the initial charge of a single Q-ball, which is expected to 
be of the same order for positive and negative charge Q-balls, and Uq {fiq) is the number 
density of positive (negative) charge Q-balls. Here, the negative charge Q-balls only appear 
in the case of e ^ e c . Notice that almost all the baryon asymmetry is initially stored in the 
Q-balls m . 



First we assume that the Q-ball decay does not produce a large extra entropy. As can be 
seen in Eqs. ( |T6| ) and (|ITD, the AD mechanism with non-renormalizable superpotentials can 
naturally provide the empirical baryon asymmetry without any extra entropy production, if 
the reheating temperature of the inflation is relatively low (e.g. Tr ~ 100 GeV for the n = 4 
case with M ~ 10 26 GeV, and for the n = 6 case with M ~ M p \.) In these cases, the number 
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density of the Q-ball is directly related to the baryon asymmetry in the present universe. It 
is easy to see that the energy density pq of the Q-ball is much less than that of the radiation 
p rad for T >T d : 

PQ _ P<t> 



Prad Prad 

_ 3 _ x m$ /ub_ 
" 4 6 TVs 



o 



~ 10 " x £ ~' (iw) (w) (^) « 1 for r >^ 1M * v - ( M > 

Therefore, the assumption of no large extra entropy production at the Q-ball decay is justified 
for e > 1CT 4 x (1 MeV/Trf). From Eq. (|45|) , we can expect that e = (9(0.1) when we use 
non-renormalizable operators in the superpotential to lift the flat direction. 
The production rate of the LSP per unit time per unit volume is given by 



dn x \ 



1 for e c < e < /3 , 



m = { €c (57) 

— for e < e c , 
e 

where N x is the number of LSPs produced per baryon number, which is at least jV™ m = 
(ef(e))- 1 > 3. 

The evolution of the number density of the LSP is then described by the following Boltz- 
mann equation: 

n B \ 6( T d-t) £ ,^ ,_..\_2 



ri x + 3Hn x = N x s[-f) v a r ' f{e) - (crv) n x , (58) 

where (crv) is the thermally averaged annihilation cross section of the LSP.Q Here, we have 
neglected the effect of the pair production of the LSPs, which is suppressed by a Boltzmann 
factor exp(— m x /T) for T < m x . 

If the Q-ball decay produces significant entropy, such as the AD mechanism without 
superpotential as discussed in Section [11 B| , uq is not directly related to the present baryon 



10 The LSPs are likely to be in kinetic equilibrium at least for T> O(MeV) pi. 
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asymmetry. In this case, the evolution of the number density of the LSPs n x is governed by 
the following coupled equations: 

ri x + 3Hn x = N x T Q n^ 1 - (av) n\ , (59) 

where ng tal is given by 

n Q total + 3Hn% ul = (for t<r d ), 

to 

l Q 

and the Hubble parameter H is obtained from the Friedman equation 



< tal = (for t>r d ), (60) 



H 2 = T^T (PQ + PX + Prad) , (61) 



where 



PQ=P<t> = {ef{e)Y l m^ (Q { - Y Q t) ng taI , (62) 

Px = m x n x ( 63 ) 
p rad + 4#p rad = ((e/(e))- 1 m ( ^ - N x m x ) Y Q n^ 1 + m x (av) n\ . (64) 

Although the Boltzmann equations have complicated forms, the final abundance of the 
LSPs can be approximately expressed by a simple analytical form. Note that the Boltzmann 
equations for the LSP [Eqs. (|58|) and (|59"|) 1 are reduced to the following one for t > r d : 

ri x + 3Hn x = — (av) r? x (for t > To) . (65) 

It may be useful to rewrite this equation in terms of Y x = n x / s and the temperature T: 



dY x 8iT 2 g* / T dgA 2 

W = ^{ 1 + 7^) {av)MplY *- m 

We assume that the s-wave contributions dominate the annihilation cross section of the 
neutralino, which is a reasonable approximation for H- and H^-like LSP with T d <C m x . By 
using the approximations g*(T) ~ g*(T d ) ~ const and (av) (T) ~ const, it can be solved 
analytically: 



Y X (T) 



Y x (T d ) V 45 



87r 2 g*(T d ) 

■ (av) M pl (T d - T) 



(67) 
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If initial abundance Y x (Td) is large enough, the final abundance Y x q for T <C is given 

by 



Y x0 ~ y a pp rox 



— (<7«) M pl T d 



(6* 



Therefore, in this case, the final abundance Y x0 is determined only by the Q-ball decay 
temperature T d and the annihilation cross section of the LSP (cry), independently of the 
initial value Y x {Tj) as long as Y X (T^) ^> Y^PP rox . In terms of the density parameter Q x , it is 
rewritten as 



h J V100 GeV 



where h is the present Hubble parameter in units of 100 km sec _1 Mpc _1 and il x = p x j p c . 
(p x and p c are the energy density of the LSP and the critical energy density in the present 
universe, respectively.) 

In the case of Y x (T d ) < y^pp rox ; the final abundance is given by 

Y x0 c W)>^(v) • ( 7 °) 

This is the case for the LSP whose annihilation cross section is small enough, such as the 
bino-like LSP. In this case, there is a very interesting feature that the relic abundance of the 
LSPs is directly related to the observed baryon asymmetry Unfortunately, however, the 
relics of the bino-like neutralino overdose the universe unless we assume an extremely light 
bino. This is easily seen from the following density parameter: 

nx >e-.(^)n B , (71) 

where m p is the mass of the nucleon. Therefore, the bino mass should be lighter than 

m x < 1 GeV x ( — )( -^-] . (72) 

We have numerically solved the Boltzmann equations for the case without a large entropy 
production [Eq. and for the case with a large entropy production [Eqs. (^)-(|64|)]. 

Here, we have used the following set of parameters: N x = 10, e = 0.1, (n B /s) = 0.7 x 10~ 10 , 
= 1 TeV and m x = 100 GeV for the former case, and iV x = 100, e < e c = 0.01, 
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1710 = 1 TeV and m x = 100 GeV for the latter case. As can be seen from Figs. [T|-[|, 
the numerical calculations reproduce the simple analytic estimation in Eq. fl68|) quite well. 
Therefore, we use the analytic estimation in Eqs. ( p8f) and (p9|) for the relic abundance of 
the neutralinos in the remainder of this paper, since it gives us the correct relic abundance 
of the LSPs as long as it is the required mass density as dark matter. 

Before closing this section, we should comment on the distribution of the LSPs. So far, we 
have assumed that they are uniformly distributed after they are produced by the Q-ball decay. 
One might wonder if the LSPs are localized near the Q-ball, since the Q-ball is a localized 
object. If this is the case, the pair annihilation rate of the LSP is enhanced and the final 
abundance might become smaller. Here, we show that this is not the case. As can be seen in 
Eq. (]67D , the abundance of LSPs approaches its final value only after (Tj, — T)/Td — 0(1), 
which means it takes a time scale At ~ r^. (This is independent of the local abundance, as 
long as it is large enough.) By that time, LSPs have spread out from the decaying Q-ball by 
a random walk colliding with background particles. Then the LSPs produced from a single 
Q-ball form a Gaussian distribution around that Q-ball, whose radius is given by f ~ \Jv7d-, 
where v~ l ~ G 2 F m x Tj 0. (Gf is the Fermi constant.) The number of the Q-balls within this 
radius is given by (4vr/3)f 3 n Q ~ 10 10 x (T d /1 GeV)- 6 (Q l /10 20 )- 1 (>x/ 1 00 GeV)~ 3 / 2 , which 
is much larger than 1. Therefore, the assumption of an uniform distribution is justified. 

V. DIRECT AND INDIRECT DETECTIONS OF HIGGSINO AND WINO 
NON-THERMAL DARK MATTER FROM AFFLECK DINE BARYOGENESIS 

In this section, we investigate the prospects of direct and indirect detection of non-thermal 
dark matter resulting from the late-time decays of Q-balls in several SUSY breaking models. 
The requirement that the LSPs produced via the late-time decays of Q-balls do not overdose 
the universe leads to the parameter region where the annihilation cross section of the LSP is 
substantially large. The possible LSP candidates are higgsino H and wino W.f^\ 

The most promising way to confirm the existence of neutralino dark matter is given by 
a direct detection through elastic scatterings of neutralino with matter. The interactions 



Here, we do not mean the pure higgsino and wino LSP. A significant fraction of the bino compo- 
nent is possible. In fact, this is the case for most of the parameter regions in the mSUGRA scenario 
we will present in this section. In the case of large tan/3, even the bino-like LSP is possible to 
achieve the desired mass density of dark matter. 
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of neutralino with matter are usually dominated by scalar couplings for relatively heavy 
nuclei A ^ 20 |41j42] . These interactions are mediated by the light h and heavy H Higgs 
exchanges, or a sfermion / exchange. Notice that the former interactions contain h\X an d 
Hxx couplings, which are suppressed for bino-like LSPs. In the case of the wino-like dark 
matter, these couplings are enhanced by the factor g2/(gl tan 9y/). As for the higgsino- 
like dark matter (except for the case of the pure higgsino LSP), these couplings are highly 
enhanced by the mixing with the gaugino components in the LSP. As we will see, these 
effects give us an intriguing possibility to detect these non-thermal dark matter in next 



generation direct dark matter searches, such as CDMS g3|, CRESST gj, EDELWEISS 
II fH, GENIUS |H and ZEPLIN MAX g7[. 

Many other indirect detection methods open up if a significant portion of halo dark 
matter consists of the neutralino LSP. These indirect methods utilize the fact that neutralinos 
accumulated in the halo or in massive bodies such as the Earth or the Sun may annihilate, 
resulting in ordinary particles, which can be detected. For example, fluxes of antiprotons 
and positrons, which are not produced in large quantities by cosmic rays, can be enhanced 
by the annihilation of neutralinos accumulated in the halo. Although the fluxes could be a 
measurable size, a clear discrimination from the background is difficult, because of the rather 
featureless spatial and energy distributions of these antiparticles. A much better signature is 
provided by neutralino annihilation into neutrinos near the centre of the Sun or the Earth. 
Since the ordinary solar neutrinos only have energies at most of the order of MeV, a multi- 
GeV neutrino signal from the Sun and the Earth may give us unmistakable evidence, although 
there are always isotropic and anisotropic backgrounds coming from neutrinos created on the 
other side of the Earth and in the outer region of the Sun by cosmic rays, respectively. 

In the present work, we investigate another promising way of detecting neutralinos in the 
halo. An excellent signature is provided by neutralino annihilation into the 27 final state 
through one-loop diagrams [[TTJ-Q This annihilation channel leads to monoenergetic 7-rays 
with energy ~ m x . Although there is an extragalactic 7 ray background and backgrounds 
of gamma-like hadronic and electron showers, the signal will stand out against them in 
favourable circumstances. In fact, this is the case for the present scenario. As we will see, a 
large fraction of H or W component in the LSP required from the late-time Q-ball decays 



12 Neutralino annihilation into the Z7 final state has a similar signature. For a detailed discussion 
of the indirect detections utilizing the neutralino annihilation channel including a photon final state, 
see Ref. g8|. 
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significantly enhances the annihilation cross section via diagrams including a chargino-W 
boson loop, which gives us an intriguing possibility to observe the 7-ray lines in the next 



generation of air Cherenkov telescopes observing the galactic centre, such as VERITAS [f|9 
HESS H and MAGIC gl] . 



In this section, we specify the parameter region where the neutralino LSPs, produced non- 
thermally through late-time decays of Q-balls, give a desired mass density of dark matter. 
In the estimation of the annihilation cross section of the neutralino, we have used all the 



tree-level annihilation channels with non-zero s-wave contributions |52],|42|]. Here, we have 
neglected the possible co-annihilation effects with the lightest chargino, which is justified as 
long as the decay temperature of Q-balls is lower than the mass difference between the LSP 
and the lightest chargino. As we will see, this condition is satisfied in most of the parameter 
region. Anyway, we are not interested in the difference by a factor of 0(1) in the estimation 
of the relic abundance, since there are ambiguities in the estimation of the decay temperature 
of Q-balls. 

We also perform calculations of the neutralino-proton scalar cross section op and anni- 
hilation cross section into the 27 final states o"2 7 in several SUSY breaking models: these 
include the "focus point" |I| in a mSUGRA scenario, the anomaly-mediated SUSY break- 



ing model || with additional universal soft scalar masses, and the no-scale model with non- 
universal gaugino masses |§ . Finally, we comment on the possibility that the non-thermally 
produced LSPs via Q-ball decays form the warm dark matter. 

A. Parameter space and possibility of detection in a mSUGRA scenario 

In the framework of minimal supergravity (mSUGRA), there are four continuous free 
parameters and one binary choice: 

m , Mi/2, Aq, tan/3, sign(/i) , (73) 

where mo, Mi/ 2 , A are the universal soft scalar mass, gaugino mass, and trilinear scalar 
coupling given at the GUT scale M G ~ 2 x 10 16 GeV, respectively. All the couplings and 
mass parameters at the weak scale are obtained through the renormalization group (RG) 
evolution. In our work, we use the SOFTSUSY code []53| to calculate quantities at the 
weak scale, which include two-loop RG equations, one-loop self-energies for all the particles 
and one-loop threshold corrections from SUSY particles to the gauge and Yukawa coupling 
constants following the methods of Ref. |53 . 



25 



At the weak scale, the higgsino mass parameter is determined by the condition of 
electroweak symmetry breaking, which, at tree level, is given by 



1 '"//„ ~ '//7/."" ' 2 >' ,,2. .,_> .,2 

2 mz ~ tan 2 /3-l 



2 nj tl u " 2 2 2 fnA\ 

m z = T 2a M — -m Hn ~ V > ( 74 ) 



where the last relation holds for moderate and large tan/3 ^ 5. 

Unfortunately, in a mSUGRA scenario, the bino-like LSP is realized in most of the pa- 
rameter space where even the thermal relics of the bino LSPs overdose the universe. In such 
a region, of course, the bino LSPs non-thermally produced through decays of Q-balls only 
make the problem worse. 

However, there exist an interesting region where the if-like LSP is naturally realized. The 
crucial observation is the existence of a "focus point" behaviour in the weak scale value of 



the soft scalar mass of the up-type Higgs multiplet, mj^ JIQ"| . The weak scale value of 



remains the weak scale even for multi-TeV mo, as long as Mi/ 2 and A are set to be around 
the weak scale. Therefore, from the relation in Eq. (|74|), one can see that the electroweak 
scale is insensitive to the input parameters for relatively large tan/3. By virtue of this focus 
point behaviour, multi-TeV values of mo do not require strong fine tunings on the input 
parameters and are natural in that sense. Such large values of m give positive contributions 
to m# , making it less negative. This, in turn, leads to smaller values of \y\ at the weak scale. 
Therefore, as mo increases further, the H content in the lightest neutralino x increases, until 
it finally enters the region \fi\ ^ 105 GeV, which is excluded by chargino searches at LEP 
II |5J. 

In the focus point region, the large H content of the LSP enhances the neutralino annihila- 
tion cross section into the W bosons via chargino exchange, which is not helicity-suppressed. 
This naturally results in a desired mass density of dark matter via the late-time decays of 
Q-balls. 

In Fig. |5], we show the allowed region for tan/3 = 15 in the (rao-Mi/2) plane. Here we 
take the sign of fi to be positive and Aq = 0. As for the criterion to select the region, we 
have imposed the following conditions: 

0.05 < tt x h 2 < 0.5 , 
1 MeV < T d < 10 GeV . (75) 

We will use these conditions throughout this paper. 

In the red shaded region, the non-thermally produced LSPs result in a cosmologically 
interesting mass density. In the allowed region, the content of the LSP is dominated by 
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H. A typical decay temperature of Q-balls, which leads to the desired mass density of dark 
matter, is 100 MeV ^ Td ^5 (a few) GeV. On the other hand, we have confirmed that the mass 
difference between the if-like LSP and the lightest chargino is at least O(10) GeV. Hence, 
the co-annihilation effects are safely neglected. This is also true for the case of tan/? = 40, 
which will be discussed later. The electroweak symmetry breaking does not occur in the 
black shaded region. The region below the blue (thick) line is excluded by the chargino mass 
limit, m x ± ^ 105 GeV. The black (thin) lines are the contours of the light Higgs boson mass, 
which are 117, 120 and 122 GeV, respectively. We also calculate the branching ratio of the 



b — ► 57 transition [56]. We adopt the following constraints on the b — > 57 branching ratio 



obtained by the CLEO experiment 

2 x 10~ 4 < B(B -> X sl ) < 4 x 10~ 4 . (76) 

In Fig. |], there is no region excluded by these bounds. 

In general, the matter-^ cross sections are dominated by spin-independent couplings for 
relatively heavy nuclei, 20. Assuming this is the case, we can obtain the detection rate 
of the neutralinos for each detector material by scaling the proton-x cross section. In Fig. ^|, 
we show this cross section in pb units in a mSUGRA scenario with tan/3 = 15. Each point 
in this scattered plot corresponds to one parameter set in the red shaded region above the 
chargino mass bound in Fig. |^. In this calculation, we adopt the following values of the 
proton matrix elements for each of the three light quarks: 

f Tn = 0.019, f Td = 0.041, f Ts = 0.U, (77) 

where fc g = (p|m g gg|p) /m p . For details about the calculation of the proton-^ cross sec- 
tion, see Refs. [ jSSp^L One can see from the figure that the proton-x cross section satisfies 
cp ~ 10~ 8 pb in most of the parameter space, which is within the reach of various next- 
generation detectors [43|-[l7]]. 



Finally, we show the annihilation rate of neutralinos into the 27 final states in Fig. [7]. Here 
we use the result of a full one-loop calculation presented in Ref . O] . The large H component 



of the neutralino enhances the annihilation rate, which even reaches 2v o<i~ f ~ 10 28 cm 3 sec 1 
in regions where the LSP is almost pure H. The fact that the if -like neutralino gives such a 



large annihilation rate was already known and investigated in Refs. [TT,5^[. However, in the 
previous works, the thermal relic of neutralinos was assumed to provide an appropriate mass 
density of dark matter. Therefore, the if-like neutralino only appears in the region with very 



large mass, m x ^ 500 GeV, because of a large annihilation cross section. In contrast, the 
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late-time decays of Q-balls lead to relatively light neutralino dark matter with a significant 
H component. This is much preferable regarding naturalness as well as detection. 

The 7-ray flux in a given direction of the sky is obtained by integrating the contributions 
along the line of sight (l.o.s.). The result highly depends on models of the density profile of 
the halo. Fortunately, all the model dependences can be factored out in terms of the following 
dimensionless function: 



8.5 kpc V0.3 GeV cm A ) Ao.s. 

where ip is the angle between the direction of the galactic centre and that of observation; 
p(l) denotes the density profile of the halo along the l.o.s.. The 7-ray flux is written by this 
function, the neutralino mass and the cross section: 

d> 7 (V>) * 1.87 x 10" 11 L/ 7 ^ A (™G?L) 2 . m cm"* sec- 1 sr" 1 , (79) 

where iV 7 is the number of photons in the final state, and iV 7 = 2 in the XX ~^ 77 annihilation. 

Although the maximum flux will be obtained in the direction of the galactic centre if) = 0, 
the experimentally relevant value is the integral of J(ip) over the solid angle around ip = 
determined by the angular acceptance of a detector. Therefore, the relevant function we have 
to treat is 



(j(o))(An) = — j^j(^)dn, (so) 

where AQ is the angular acceptance of the detector. The value of ( J(0)) and its dependence 
on AQ in various halo models are discussed in Ref. |59[| . A typical value is AQ = 10~ 3 sr 
for a generic next-generation air Cherenkov telescope (ACT). We define the value of (J{0)) 
averaged over this angular acceptance as j = («/(0)) (10~ 3 sr). 

In Fig. [?|, we also plot 5a sensitivity curves for typical next generation ACT arrays |49]-|5T 



observing the galactic centre, adopting several models of dark matter distribution, including 
a Moore et al. profile |60| for j = 10 5 and Navarro, Frenk and White profile |H] for j = 10 3 . 



In obtaining the sensitivities curves, we consider an instrument with a 10 9 cm 2 effective area 
and 15% energy resolution, assuming 10 6 sec exposure and using the standard estimation 
of the background presented in Ref. We have also adopted a hadronic rejection that 



is improved by a factor of 16 with respect to the current Whipple detector, following the 
arguments presented in this reference. One can see that by assuming somewhat cuspy dark 
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matter profiles of the halo, we have a relatively large possibility for detecting monoenergetic 
7-rays, especially for the region with almost pure higgsino LSP. 

In Figs. D, P and 10, we show, respectively, the allowed parameter space, the proton-^ 



cross sections and ^-annihilation rates into 27 final state for the mSUGRA scenario with 
tan/5 = 40. We take sign(/x) to be positive also in this case, which is desirable to avoid a 
large branching ratio of the b — > 57 transition. Conventions are the same as in Figs. |5], ^| and 
0. One can see some differences in Fig. |8] compared with Fig. |5|. First, the region below the 
green (dot-dashed) line is excluded by the constraints on the branching ratio of b — > 57 from 
the CLEO experiment, B(B — ► X s 7) > 2 x 10 -4 . Second, the allowed region is significantly 
widened and even the region with the bino-like LSP (the region for mo ^ 1 TeV) can lead 
to the desired mass density of dark matter. This is because, in the case of large tan/5, the 
annihilation cross section through the process XX ~ > A ~ y ff is strongly enhanced. The 
reason is that the coupling with fermions Aff is proportional to tan/3 and that the mass of 
the pseudo-Higgs boson is significantly reduced by large down-type Yukawa couplings, 
which results in a much smaller suppression factor (m^/m^ 4 . 

Here, we should add a comment. One may think that the formation of Q-balls does not 
occur in the focus point region, since a large soft mass squared Wq would lead to the ^-factor 
in Eq. (|33j) being positive. However, even if = O(10 TeV), the i^-factor is negative as 
long as we adopt flat directions, which are solely constructed by the first and second families. 
This condition can easily be satisfied in most of the flat directions. 

B. Parameter space and possibility of detection in the anomaly- mediation model 

In pure anomaly-mediated SUSY breaking, soft terms are determined by RG-invariant 
expressions involving the gauge and Yukawa couplings; hence the soft terms are completely 
fixed by the low energy values of these couplings and an overall scale 1713/2'- 

Mx = -YlMil 771 ^ = 7 m3/2 ' 

2 1 d 2 lnZ<s, 2 1 / <9 7 <9 7 




where g, y are gauge and Yukawa couplings, respectively, <3> denotes the general superfield, 
and m%, Zq> are its soft scalar mass squared and wave function renormalization factor; M\ 
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is the soft mass for the gaugino A, A y is the trilinear scalar coupling associated with the 
Yukawa interaction with a coupling constant y. Here the sum J2i runs through the fields 
included in this Yukawa interaction. We have defined the renormalization group functions as 
j(g,y) = dlnZ / dlia.fi, /3 g (g,y) = dg/d\n.{i, and f3 y = dy/dlrifi. 

As you can see from these relations, the soft masses squared for sleptons are negative, 
which is the biggest difficulty in pure AMSB models. Although many possible solutions are 
proposed |62|], we adopt, in the present work, a phenomenological solution to the negative 
slepton mass problem. We assume the additional universal scalar mass mo at the GUT scale, 

ml = ~\ (f^ 9 + f^) m '/ 2 + m ° ' (82) 

and then evolve the RG equations to obtain the low energy spectrum. Then the entire 
parameter space is specified by the following 4 parameters in this minimal framework: 

m 3 / 2 , mo, tan/5, sign(/i) . (83) 

From Eq. fl8"T|), we see that the ratios between the gaugino masses at the weak scale are 
approximately given by 

Mi : M 2 : M 3 « 2.8 : 1 : -8.3 , (84) 

where M%, M 2 and M 3 are the gaugino masses for U(l)y, SU(2)l and SU(3)c, respectively. 
Therefore, in the present scenario, the W LSP is realized in most of the parameter space. 

In Fig. |ll], we present a parameter space where the LSPs produced by late-time decays of 
Q-balls provide a desired mass density of dark matter and the proton-x scalar cross section 
in that region. Here we take tan/3 = 15, and sign(/x) negative to avoid too large a branching 
ratio of b —>■ 57. There is no region excluded from the constraints on b — > 57. The black 
shaded region is excluded by the f LSP and/or overdosing universe by the LSPs produced 
via Q-ball decays. In all of the remaining white space, the W LSP from the Q-ball decay can 
be the dominant component of the dark matter. The typical decay temperature of Q-balls 
resulting in the suitable mass density of dark matter is 10 MeV < T d < 100 MeV. On the 
other hand, the mass splitting between the charged and neutral Ws is of order 100 MeV to 
1 GeV 0,0, which allows us to neglect the co-annihilation effects. 



We present contours of the proton-x cross section by the red (thick) lines, which corre- 
spond to ap = 10~ 8 , 10~ 9 , . . . , 10~ 12 pb from left to right, respectively. A strange struc- 
ture appearing in the lower right-hand corner indicates a contamination of the H compo- 
nent in the LSP. The green (dotted) lines are the contours of the lightest neutralino mass, 
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m x = 100, 150, . . . , 300 GeV, respectively. We also present the contours of heavy Higgs 
boson mass by the blue (dashed) lines, which denote ms = 500, 750, . . . , 2000 GeV from 
left to right, respectively. In contrast to the previous expectation ^,65], the direct detection 
rates are rather small, and only a restricted region is within the reach of the next-generation 
detectors, o> ^ 10~( 9-8 ) pb. This comes from the large mass of the heavy Higgs boson rxiu- In 
relatively large tan/5, the proton-x cross section is dominated by H exchange, and it reduces 
as o P oc mjj 4 . In Refs. |2|,[55||, it was assumed that m H ~ 300 GeV, which is realized only in 



a small space, as one can see from Fig. [11 



Fortunately, we may obtain a clear signal from the neutralino annihilation into the 27 
final state [pq] . In Fig. 12, we show the annihilation rate of the LSP into the 27 final state. 
Each orange dot falls in somewhere in the allowed space presented in Fig. 11. The blue lines 
are 5a sensitivity curves for j = 10 2 and 10 3 , which are estimated by the same parameter 
set as was used in Figs. |7| and pi], but we do not assume any improvements on the hadronic 
rejection and use ehad = 1 in Ref . ]5]| . Even if we do not assume strongly cuspy profiles of the 
dark matter density in the halo, a large portion of the parameter space predicts a measurable 
size of the 7-ray flux. 

In Figs. [13] and [14|, we show the plots for the case of tan/3 = 30 and sign(/i) < 0. 



Conventions are almost the same as those in Figs. [TT| and [12], respectively, but they are 
explicitly denoted in the caption to each figure. In such a large tan/3, we can see that 
the proton-;^ cross section is strongly enhanced. This is because the if-nucleon-nucleon 
coupling, which dominates the proton-x cross section, is proportional to tan/3, and the large 
down- type Yukawa couplings reduce the mass of the heavy Higgs boson m 2 H . We can expect 
the direct detection of the W dark matter in a relatively wide parameter space in the next 
generation of detectors. In Fig. [y| a small region with m ^ 500 GeV, m 3 / 2 ~ 30 TeV is 
excluded by the constraint B(B — > X s j) > 2 x 10~ 4 p7| , although it is not depicted in the 
figure. The neutralino annihilation rate into the 27 final state is almost the same as in the 
case of tan/? = 15, which is consistent with the result of Ref. Q . 

In AMSB models, the neutralino annihilation into the Z'j final state is known to be nearly 
twice that into the 27 final state (see Ref. Note that, in such a minimal framework of 

AMSB models, the thermal relic of the W cannot be a significant component of dark matter, 
because of its large annihilation cross section. Such a large annihilation cross section of the 
LSP is, on the contrary, much more advantageous for the AD baryogenesis with late-time 
decays of Q-balls, although the conditions given in Eq. ([IT]) or (^) should be satisfied. 
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C. Parameter space and possibility of detection in no-scale models 
with non-universal gaugino masses 



In no-scale models [p8| , all the soft parameters except gaugino masses are assumed to 
vanish at some high energy scale Mj, which is usually taken to be the GUT scale Mx = 
Mq- The soft parameters at the weak scale are generated by RG effects dominated by 
the gaugino masses, which are automatically flavour-blind and naturally suppress the SUSY 
flavour-changing neutral current (FCNC) interactions Recently, models with the no-scale 
boundary condition have attracted much attention, since a natural realization was found as 
gaugino-mediation models [7C . 

Unfortunately, the MSSM with the no-scale boundary condition at the GUT scale was 
found to be inconsistent with phenomenological requirements p,|7l,72|. This is mainly due 
to the lower bound on the light Higgs boson mass and the cosmological requirement that 
a charged particle must not be the LSP. Several solutions were proposed to this problem, 
which are, imposing the no-scale boundary conditions above the GUT scale [ff3| , gauging the 
U(l) b-l symmetry |73| and assuming the non-universal gaugino masses at the GUT scale 0. 

In the present work, we pick up the last scenario with non-universal gaugino masses. In 
this scenario, as investigated in Ref. ||, the W LSP is realized in a wide parameter space, 
which is expected to derive a desired mass density of dark matter from late-time decays of 
Q-balls. 

In Fig. [TjJ we present allowed regions where the resultant LSPs give a desired mass density 
of dark matter and the proton-x cross section in that region. Here, we take M 2 = 200 GeV, 
tan/3 = 10 and vary the ratios of the gaugino masses at the GUT scale Mi/M 2 and M 3 /M 2 , 
which are denoted by the x-axis and the y-axis, respectively. The red (thick) solid lines are 
the contours of the proton-;\; cross section o>, which are Up = 10 -7 , 10~ 8 and 10 -9 pb from 
the bottom up. The black shaded region is excluded, since the f is the LSP, or the resultant 
LSPs overdose the universe, or the EWSB cannot be implemented. The blue (dashed) lines 
are the contours of the Higgs boson mass, nth = 114.1 and 120 GeV, and the region below 



the lower line is excluded by the Higgs boson mass bound [75|. The region below the green 
(dot-dashed) line is excluded by the constraints on the b — > 57, B{B — > X s ^) > 2 x 10~ 4 . 
In Fig. [16], we show the neutralino annihilation cross section into the 27 final state for the 



allowed region in Fig. 15. We also present the 5<r sensitivity curve for j = 10 . Conventions 



are the same as those in Fig. [L2| for anomaly- mediation models. In the allowed region, the LSP 
is given by the PU-like neutralino, which predicts large annihilation cross sections comparable 
with those in anomaly-mediation models. 
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In Figs. and [Tj|, we present the corresponding figures in the case of tan/3 = 30. Con- 
ventions are the same as those in Figs. |l^ and [16|. (M2 = 200 GeV is assumed at the GUT 
scale.) In this case, the direct detection rate is enhanced by large tan/3, but a large por- 
tion of the parameter space is excluded by the constraint from the b — > 57. The bino-like 
LSP is realized in the small spot appearing on the left of the excluded region. In this re- 
gion, the large annihilation cross section required from the late-time decays of Q-balls is 
provided by the pseudo-Higgs A exchange diagram via the small H contamination in the 
LSP, whose amplitude is enhanced by the large tan/3 = 30 and the relatively small tua- 
In this bino-like LSP region, the annihilation rate into the 27 final state is relatively small 
2va 2l ~ 3 x 10~ 29 cm 3 sec -1 . 

Warm dark matter production from late-time decays of Q-balls 

Before closing this section, we should remark on an interesting possibility that the resul- 
tant H-like LSPs (or maybe 5-like LSP in the case of a large tan/?) form the warm dark 
matter (WDM). As we have seen, cuspy structures in the halo profiles have a drastic conse- 
quence on the future observations in indirect dark matter searches. However, it is argued that 
the neutralino cuspy profile might be inconsistent with the radio emission from the centre of 
the Galaxy [J7fj , although further detailed studies are needed to fix the situation. Maybe the 



inflation models with a suppressed strength in the small-scale perturbations might solve the 
problem. Another possible solution is the WDM. 

If the neutralino LSPs are produced after the decoupling of thermal interactions with rela- 
tivistic velocity, they may serve as a WDM and wash out the cuspy profiles in the halo [|77|,[78] • 



The conditions needed for the MSSM neutralino to form a WDM were studied in Ref. P0 
Actually, these authors found that it is possible if the temperature is low enough when the 
neutralinos are produced via decays of non-thermal sources: 

T d < 5 MeV for E-like LSP 

and 

T d <2MeV for H-\ike LSP. (85) 

As for the W LSP, they found it is impossible to form the WDM. If the above condition is 
satisfied, energy reductions of neutralinos by scatterings with thermal background are small 
enough. Such a situation may arise in some models, for example, the late-time decays of 



33 



the heavy moduli 0,^] and the evaporation of cosmological defects [78|. Here, we propose 



another way to generate WDM via the late-time decays of Q-balls. 

To obtain the sufficient free streaming length needed to suppress the small-scale pertur- 
bations, the current velocity of the LSP should satisfy [ff8| : 



v Q ~ 1(T (7 ~ 8) . (86) 

If neutralinos are produced via late-time decays of Q-balls, the current velocity of the neu- 
tralinos is given by 

Vo = ^ K 2 . 4 x 10 -» (vm (-^\ (m^) , m 

T d m x \ T d / V10 TeV/ y m x J 

where T = 2.7 K is the current temperature of the cosmic microwave background and T d is 
the decay temperature of Q-balls; is the mass of the relevant flat direction field <fi. Here, 
we have assumed that there is no substantial energy reduction by scattering with thermal 



backgrounds, which is justified if the condition in Eq. (85) is satisfied. 

Therefore the resultant neutralinos from late-time decays of Q-balls may serve as a WDM, 
for example, in "focus point supersymmetry" with m = 0(10 TeV) or in "effective SUSY 
models" with very large soft masses for the squarks in the first and second families m~ = 
0(10 TeV) 2 |8|. 

VI. CONCLUSIONS AND DISCUSSIONS 

Both of the production mechanisms of the observed baryon asymmetry and natures of 
dark matter are among the most fundamental problems in particle physics as well as in cos- 
mology. In the SUSY framework, there exist rather stringent upper bounds on reheating 
temperatures of inflation because of the "cosmological gravitino problem", which strongly 
constrains various baryo/leptogenesis scenarios requiring relatively high reheating temper- 
atures. Affleck-Dine baryogenesis is one of the most promising candidate to generate the 
observed baryon asymmetry in low reheating temperatures, which is free from the gravitino 
problem. 

In recent developments, it has become clear that the formation of Q-balls and their 
late-time decays are almost inevitable consequences of AD baryogenesis.^] This requires a 



13 



This is not the case for the AD leptogenesis via the LH U flat direction. See Sec. II. 
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significantly large annihilation cross section of the LSP for the resultant LSPs not to overdose 
the universe. 

In this paper, we investigated in detail the predictions of the late-time decays of Q-balls. 
First, we reviewed the AD baryogenesis, the Q-ball formation and decays, with thermal effects 
taken into account. Second, we discussed the non-thermal production of dark matter via the 
late-time decays of Q-balls with and without entropy production. Finally, we specified the 
allowed parameter space where the neutralinos produced via the decays of Q-balls result in a 
desired mass density of the dark matter in several SUSY breaking models, and then discussed 
the prospects of direct and indirect detection of these non-thermal dark matter candidates. 
We also discussed the possibility of WDM generation via Q-ball decays. 

As an interesting new twist, we pointed out that the AD baryogenesis without superpo- 
tential is now one of the most beautiful scenarios. Quite a reasonable set of parameters can 
explain the observed baryon asymmetry without any assumptions on additional late-time 
entropy productions, such as the decay of heavy moduli or a thermal inflation. The resultant 
baryon asymmetry and dark matter density are determined solely by the potential for the 
AD field and the annihilation cross section of the LSP. They are completely independent of 
the reheating temperature of the inflation and any other details in the history of the universe. 

The parameter space and detection possibility for each SUSY breaking model can be 
summarized as follows: 

• In minimal supergravity scenarios, allowed parameter space appears in the "focus-point" 
region where the dominant component of the LSP is provided by H. A large content of H in 
the LSP strongly enhances the direct detection rate through the Higgs exchange diagrams. 
As a result, the direct detection rate is within the reach of various next-generation detectors 
in most of the parameter space. The detection rate of the monoenergetic 7-rays via neu- 
tralino annihilation is also enhanced with respect to bino-like dark matter, especially in pure 
H LSP regions. Somewhat cuspy density profiles of the halo j ^ 10 35 allow us to observe the 
monoenergetic 7-ray lines in the case of pure H dark matter. 

• In anomaly-mediated SUSY breaking models, the LSP is VU-like in most of the param- 
eter space, which is now a promising candidate of the dark matter because of the late-time 
decays of Q-balls. The direct detection rate highly depends on tan/5 and becomes very large 
in the case of large tan/3. This is because the coupling of the heavy Higgs boson with the nu- 
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cleon is proportional to tan/3, and also the mass of the heavy Higgs boson is strongly reduced 
by the large down-type Yukawa couplings. The indirect detection rate using monoenergetic 
7-rays is surprisingly large, as reported in Ref. HB"6"| , and is almost independent of tan/3. Even 



if we do not assume a highly cuspy density profile, the detection rate is within the reach of 
the next generation of air Cherenkov telescopes [fj5H °Tl • 



• In no-scale models with non-universal gaugino masses, the LSP is HMike in a wide 
parameter space. (The region with the H LSP is excluded by the Higgs mass bound and 
the constraints on the h — > sj branching ratio.) The direct and indirect detection rates are 
comparable with those in anomaly-mediation models, and within the reach of various next- 
generation detectors. 

Although we have concentrated on the monoenergetic 7-rays for the indirect detection in 
this paper, the possibility of detecting the neutralino in other indirect dark matter searches 
is also significantly enhanced in the case of H and W dark matter. The prospects of other 
indirect detections in mSUGRA scenarios can be obtained in Refs. ||81|| . As for the anomaly- 



mediation models, see Refs. |66|j65|] . Last year, the recent HEAT experiment |82j has con- 



firmed an excess of a high energy positron flux in cosmic rays ||33| . If the LSP is H- or 



IV-like, the annihilation into a pair of W-bosons always dominates as long as m x > mw 
The subsequent decay of a W-boson into a positron, W + — > e + + u, will produce a large 
excess of positron flux at and below an energy of about half the W-boson mass. In fact, it 
was shown in Ref. pi that the observed positron spectrum is naturally explained by the 



annihilation of neutralino dark matter, if it is H- or IV-like non-thermal dark matter. This 
fact may indicate the existence of non-thermal sources for these LSPs in the history of the 
universe. 

Finally, we comment on some advantages of our scenario relative to the others that can 
produce the H or VT-like dark matter. Actually, there exist several other mechanisms that 
can generate the H- or PU-like non-thermal dark matter, which are the decays of heavy moduli 
fields and/or heavy gravitinos with a mass of the order of 100 TeV J55| |, and the evaporations 



of topological defects []7q]. The former case only appears in the anomaly- mediated SUSY 
breaking scenarios. Furthermore, the decays of heavy moduli fields are usually accompanied 
by an extra large entropy production, which causes a serious problem to explain the observed 
baryon asymmetry. If one relies on the AD baryogenesis to produce enough baryon asymme- 
try, then the resultant late-time Q-ball decays also play a role as the non-thermal source for 
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the lU-like dark matter. As for the evaporation of topological defects, we have no clear reason 
to expect that they take place at the appropriate temperature of 1 MeV^S Td ~( a f ew ) GeV, 
and that they produce the enough initial abundance of the H- or the lU-like LSP to be a 
dominant component of dark matter. They may also cause a serious problem by diluting the 
baryon asymmetry. 

In our scenario, the reasonable set of parameters in the scalar potential of the AD field 
to explain the observed baryon asymmetry naturally leads to the desirable decay tempera- 
ture of the Q-balls for the non-thermal production of these LSPs in various SUSY breaking 
scenarios. The baryon asymmetry is, of course, guaranteed by the AD mechanism. 

As we have seen in this article, H- and lU-like LSPs are now promising candidates for the 
dark matter. The possibility of direct and indirect detection of these dark matter candidates 
is much larger than that of the standard bino-like LSP. If H or W dark matter is indeed 
confirmed in future experiments, it may shed a bright light on the origin of the whole matter 
in our universe. 

Finally, we add a brief comment on the fine tunings on the soft SUSY breaking parameters 
to explain the mass density of dark matter for a given decay temperature of a Q-ball. As 
we have discussed in the text, the dominant annihilation channel of the higgsino- and wino- 
like LSP is given by the decay channel into a pair of gauge bosons via the lightest chargino 
exchange. This means that the annihilation cross section of these LSPs are rather stable 
under the variation of other soft parameters, such as slepton and squark masses. (Note that 
the variation of the squark masses has only a mild effect on the size of /x-term in the focus 
point region.) Consequently, no strong fine tuning on the soft SUSY breaking parameters is 
required for a given decay temperature of a Q-ball. This is a clear contrast to the standard 
pure bino dark matter, which needs a strong degeneracy between the bino and the lightest 
stau. 
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FIGURES 
T d = 10 MeV (without entropy production) 
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FIG. 1. The evolution of the abundance of the neutralino dark matter generated from the Q-ball 
decay for Td = 10 MeV with (av) = 10 -7 GeV~ 2 and 10 -5 GeV -2 , which are represented by thick 
solid lines. The abundances estimated by the analytic formula in Eq. ( |68|) are shown in dashed lines. 
In this figure, we have assumed that the energy density of the Q-ball is small enough with respect 
to that of the radiation. The parameters are taken to be = 1 TeV, m x = 100 GeV, e = 0.1 and 
iV v = 10. 
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FIG. 2. The same as Fig. [T], but with T d = 1 GeV and (av) = 10" 
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T d = 10 MeV (Q-ball dominant) 

e<e=0.01, N=100, m=lTeV, m =100GeV 
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FIG. 3. The evolution of the abundance of the neutralino dark matter generated from the Q-ball 
decay for Td = 10 MeV with (av) = 10~ 7 GeV -2 and 10~ 5 GeV -2 , which are represented by thick 
solid lines. The abundances estimated by the analytic formula in Eq. ( |68| ) are shown in dashed 
lines. In this figure, we have assumed that the Q-balls dominate the energy density of the universe 
before their decay. The parameters are taken to be = 1 TeV, m x = 100 GeV, e < e c = 0.01 and 
N Y = 100. 
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FIG. 4. The same as Fig. g but with T d = l GeV and (av) = 10 -9 , IO" 7 , and IO" 5 GeV" 2 . 
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FIG. 5. The allowed region in the mSUGRA scenario with tan/? = 15 and Aq = in the 
{rriQ-Mi/2) plane. In the red shaded region, non-thermally produced LSPs via decays of Q-balls 
result in a cosmologically interesting mass density. The black shaded region is where the electroweak 
symmetry breaking cannot be implemented. The region below the blue (thick) line is excluded by 
the chargino mass bound "m x ± ~ 105 GeV. The contours of the light Higgs boson mass are given 
by the black (thin) lines, which correspond to mh = 117, 120, 122 GeV, respectively. There is no 
region excluded by the bounds on the b — > S7 branching ratio. 



44 




0.0001 



100 150 200 250 

m x [GeV] 



300 



FIG. 7. Annihilation rate of neutralinos into the 27 final state in the mSUGRA scenario for the 
red region in Fig. |5l also with the chargino mass bound. 
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FIG. 8. The allowed region in the mSUGRA scenario with tan/3 = 40 and Aq = in (mo-M^) 
plane. In the red shaded region, non-thermally produced LSPs via decays of Q-balls result in a 
cosmologically interesting mass density. The black shaded region is where the electroweak symmetry 
breaking cannot be implemented. The region below the blue (thick) line is excluded by the chargino 
mass bound fn x ± ^ 105 GeV. The contours of the light Higgs boson mass are given by the black 
(thin) lines, which correspond to = 114.1, 117, 120, 122 GeV, respectively. The region below 
the green (dot-dashed) line is excluded where the branching ratio of b — > 57 violates the CLEO 
bound, B(B X s j) > 2 x 10~ 4 . 
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FIG. 9. Proton— x cross section in the mSUGRA scenario for the red region in Fig. 
the chargino mass bound and the constraint from b — > 57. 
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FIG. 10. Annihilation rate of neutralinos into the 27 final state in the mSUGRA scenario for 
the red region in Fig. also with the chargino mass bound and the constraint from b — > 57. 
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FIG. 11. Allowed parameter space and proton-x cross section in the anomaly- mediation model 
with tan/3 = 15. Here, we take sign(/z) negative. In this plot, there is no region excluded by b — ► 57. 
The black shaded region is excluded by f-LSP or the produced neutralinos overdose the universe. 
The wide white region leads to a desired mass density of dark matter via decays of Q-balls. The 
red (thick) lines are contours of proton-^ cross section, which are op = 10~ 8 , 10~ 9 , 10~ 10 , 10 -11 
and 10 -12 pb, from left to right, respectively. The blue (dashed) lines are contours of the mass of 
the heavy Higgs boson mH, which are, from left to right, 500, 750, 1000, . . . 2000 GeV. The green 
(dotted) lines are contours of the lightest neutralino mass m x for 100, 150, . . . 300 GeV. 
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FIG. 12. Annihilation rate of neutralinos into the 27 final state in the AMSB model with 
tan/3 = 15. Each orange dot corresponds to one parameter set in the white region in Fig. 11. The 
two blue lines denote the ha sensitivity curves for j = 10 2 , 10 3 . 
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FIG. 13. Allowed parameter space and proton-x cross section in the anomaly-mediation model 
with tan/3 = 30. Here, we take sign(//) negative. The black shaded region is excluded by f-LSP or 
the produced neutralinos overdose the universe. The black region also includes small regions where 
we cannot obtain the convergence of SOFTSUSY code. The wide white region leads to a desired 
mass density of dark matter via decays of Q-balls. The red (thick) lines are contours of proton-^ 
cross section, which are ap = 10 -7 , 10~ 8 , 10~ 9 , 10~ 10 and 10 -11 pb, from left to right, respectively. 
The blue (dashed) lines are contours of the mass of the heavy Higgs boson mn, which are, from 
left to right, 500, 750, 1000, . . . 1750 GeV. The green (dotted) lines are contours of the lightest 

150, 200 GeV. 
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FIG. 14. Annihilation rate of neutralinos into the 2j final state in the AMSB model with 
tan/3 = 30. Conventions are the same as those in Fig 12. 
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FIG. 15. Allowed region and the proton-^ scalar cross section ap in the no-scale model with 
non-universal gaugino masses with tan/3 = 10. Here, we take M<i = 200 GeV at the GUT scale. The 
red (thick) lines are the contours of the up = 10 -7 , 10~~ 8 , 10 -9 pb from the bottom up, respectively. 
The black shaded region is excluded by the fact that f is the LSP, or the resultant LSPs from 
late-time decays of Q-balls overdose the universe, or the EWSB cannot be implemented. The blue 
(dashed) lines are the contours of the Higgs boson mass = 114.1, 120 GeV and the region 
below the lower line is excluded. The region below the green (dot-dashed) line is excluded since 
B(B -► X sl ) < 2 x 10" 4 . 
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FIG. 16. Annihilation rate of neutralinos into the 27 final state for the allowed region in Fig. 15 
Conventions are the same as those in Fig. [l^. 
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FIG. 17. Allowed region and the proton-^ scalar cross section dp in the no-scale model with 
non-universal gaugino masses with tan/3 = 30. Conventions are the same as those in Fig. |l5| The 
bino-like LSP is realized in the small spot appearing on the left of the excluded region. The required 
large annihilation cross section is obtained by the A exchange diagram enhanced by the large tan/? 
and the relatively small rriA- 
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FIG. 18. Annihilation rate of neutralinos into the 27 final state for the allowed region in Fig. ITj 



Conventions are the same as those in Fig |12j. The dots with smaller annihilation rates correspond 
to the region with the bino-like LSP. 
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